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Summary 


The  following  objectives  have  been  achieved: 


•  studied  the  reaction  kinetics  of  Mg+  ions  with  O3,  O2,  N2,  CO2,  H2O  and  N2O  using  a 
laser  flash  photolysis-laser  induced  fluorescence  technique 

•  studied  the  reactions  of  MgO+  and  Mg02+  with  atomic  O,  and  a  series  of  ligand-switching 
and  association  reactions  of  Mg+  cluster  ions  with  O2,  H20  and  CO2,  using  a  flow  tube- 
mass  spectrometry  technique 

•  studied  the  reactions  of  MgO,  Mg02,  Mg03  and  MgC03  with  atomic  O  using  a  flow 
tube-  laser  induced  fluorescence  technique 

•  produced  a  set  of  high-level  quantum  theory  calculations  of  the  thermochemistries  of  all 
these  reactions 

•  completed  an  analysis  of  satellite  observations  (using  the  SCIAMACHY  instrument  on 
ENVISAT)  of  Mg  and  Mg+;  this  has  shown  that  unlike  other  metals  in  the  mesosphere, 
atomic  Mg  is  not  removed  in  the  presence  of  noctilucent  ice  clouds 

•  developed  a  new  model  of  magnesium  chemistry  in  the  mesosphere/lower  thermosphere 
which  predicts  absolute  densities  of  the  Mg  and  Mg+  layers,  including  their  seasonal 
variations,  in  excellent  agreement  with  satellite  observations. 

•  applied  this  model  to  determine  the  lifetimes  of  Mg+  ions  in  sporadic  ion  layers  in  the 
terrestrial  and  Martian  atmospheres 
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magnesium  ion  chemistry  in  the  upper  atmosphere,  under  review  at  Physical  Chemistry 
Chemical  Physics. 
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Introduction 


Metallic  species  are  deposited  into  the  Earth’s  upper  atmosphere  by  the  ablation  of 
approximately  50  tonnes  of  interplanetary  dust  that  enters  the  atmosphere  each  day.1  The  dust 
undergoes  frictional  heating  to  its  melting  point  (~  1800  K),  after  which  metallic  species  ablate 
and  are  deposited  in  the  mesosphere-lower  thermosphere  (MLT)  region  of  the  atmosphere  (70  - 
120  km).2 


Magnesium  is  one  of  the  most  abundant 
metals  in  the  MLT,  with  layers  of  Mg+  being 
observed  by  rocket-borne  mass 
spectrometry3,4  and  satellite  observations  of 
the  earth’s  dayglow.5'15  Recent  satellite 
observations,  made  by  the  GOME16  and 
SCIAMACHY17,18  instruments,  have  found 
Mg  to  have  a  column  density  of  ~  2  x  109 
cm'2,  with  Mg+  exhibiting  column  densities 
between  3  x  109  cm'2  at  its  winter  minimum 
and  1.2  x  1010  cm'2  at  its  summer  maximum 
at  high  latitudes.  Mg+  has  a  peak  altitude 
between  90  and  100  km  and  a  peak  ion 
concentration  of  (1  -  5)  x  10  cm'  during 
daylight,  which  declines  to  10  cm'  at 
night.19  Mg  exhibits  lower  concentrations  and 
a  peak  altitude  around  88  km.17,20  Unlike  the 
other  prominent  meteoric  metals  (Fe,  Na,  K 
and  Ca),  neither  Mg  nor  Mg+  can  be  observed 
by  ground-based  lidar  (laser  radar)  as  they 
have  resonance  transitions  in  the  UV  region  at 
285  and  280  nm  respectively,  where  light  is 
strongly  absorbed  by  stratospheric  O3. 

Mg+  is  produced  from  Mg  by  photo¬ 
ionization  and  charge  transfer  with  NO+  and 
C>2+,  which  are  the  dominant  ions  in  the  lower 
thermosphere.  A  striking  difference  between  magnesium  and  other  meteoric  metals  is  the  large 
ratio  between  Mg+  and  Mg.  This  ratio  ranges  from  about  1.5  to  more  than  10.  Magnesium  seems 
to  be  the  only  metal  displaying  such  a  large  ion/neutral  ratio.  Na+/Na  and  Fe+/Fe  display  ratios 


Figure  1.  Schematic  diagram  of  magnesium 
chemistry  in  the  MLT.  Blue  arrows:  reactions 
studied  in  the  EOARD  project.  Black  arrows: 
previously  studied  by  the  Leeds  group.  Grey 
arrows:  remain  to  be  studied,  but  rate  constants 
can  be  estimated  reliably.  Thicker  arrows  indicate 
more  important  reaction  pathways. 
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of  ~  0.222’23  and  Ca+/Ca  has  a  ratio  of  ~  2. 24  This  is  even  more  striking  as  Mg+  is  not 
significantly  depleted  relative  to  other  metals  in  the  MLT. 

Metallic  ions,  mainly  Mg+  and  Fe+, 22,25  are  the  major  constituents  of  sporadic  E  layers  (Es).  Es 
are  thin  layers  of  concentrated  plasma  between  1  and  3  km  wide,  which  occur  at  altitudes 
between  85  and  140  km.26  They  play  an  important  role  in  over-the-horizon  and  space-to-ground 
radio  communications.  Several  mechanisms,  including  wind  shear  and  gravity  waves,  have 

96 

been  proposed  to  explain  their  formation. 

Figure  1  is  a  schematic  diagram  of  the  chemistry  of  magnesium  in  the  MLT.  The  blue  arrows 
indicate  the  reactions  studied  in  this  EOARD  project,  demonstrating  that  the  chemistry  is  now 
much  better  constrained  and  understood.  Consider  the  ion-molecule  chemistry  first.  Mg+  ions 
are  produced  directly  during  meteoric  ablation  via  hyperthermal  collisions  between  Mg  atoms 
and  air  molecules,  and  also  by  charge  transfer  between  Mg  and  the  ambient  ions  NO+  and  C>2+ 
(solar  photo-ionization  of  Mg  is  minor).  Above  110  km  the  neutralization  of  Mg+  by  dielectric 
recombination  (Mg+  +  e  —*■  Mg  +  hv )  is  the  only  significant  removal  process  of  Mg+  ions. 
However  this  process  is  very  inefficient,  with  an  estimated  rate  coefficient  of  1  x  10~12  cm3  s'1 28 
implying  an  average  lifetime  for  Mg+  in  excess  of  10  days. 

At  lower  altitudes  ion-molecule  reactions  become  important: 


Mg+  +  O3  — >  MgO+  +  O2 

AH(0  K)  =  -1 16  kJ  mol’1 

(1) 

Mg+  +  N2  +  M  — >•  Mg+.N2  +  M 

AH(0  K)  =  -30  kJ  mol'1 

(2) 

Mg+  +  O2  +  M  — *  MgC>2+  +  M 

AH(0  K)  =  -90  kJ  mol'1 

(3) 

Mg+  +  C02  +  M  — ►  Mg+.C02  +  M 

AH(0  K)  =  -67  kJ  mol'1 

(4) 

Mg+  +  H20  +  M  -►  Mg+.H20  +  M 

AH(0  K)  =  -131  kJ  mol'1 

(5) 

where  M  is  the  third  body  (principally  N2  in  the  terrestrial  atmosphere)  which  removes  excess 
energy  and  stabilizes  the  association  complex.  The  reaction  enthalpy  changes  are  determined 
from  high-level  quantum  theory  calculations  which  were  published  recently  as  part  of  this 
project.  Note  that  the  molecular  ion  nomenclature  adopted  here  is  to  indicate  an  ion  where  there 
is  a  chemical  bond  as  MgX+,  and  one  where  there  is  essentially  only  an  electrostatic  interaction 
as  Mg+.X. 

Following  the  pioneering  work  at  the  NOAA  Aeronomy  Laboratory  in  the  1960s,  formation  of 
MgO+  via  reaction  1  is  known  to  be  rapid,  and  should  dominate  Mg+  removal  above  90  km  in  the 
terrestrial  atmosphere.30,31  Below  85  km,  reaction  with  O2  should  become  more  important.30  In 
contrast,  reaction  4  dominates  in  the  CCVrich  Martian  atmosphere.32  Although  the  Mg+.N2 
complex  has  been  observed  in  the  laboratory,33  the  kinetics  of  reaction  2  do  not  appear  to  have 
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been  studied.  Reaction  5  is  comparatively  fast,34  but  H2O  is  a  trace  constituent  and  so  this 
reaction  should  not  be  competitive. 

The  molecular  ions  formed  by  these  reactions  will  eventually  undergo  very  efficient  dissociative 
recombination  (DR)  with  an  electron  (Mg+  .X  +  e — >  Mg  +  X)  to  produce  neutral  Mg.  However, 
before  this  occurs  these  ions  can  undergo  rapid  switching  reactions,  such  as: 

Mg+.C02  +  H20  -»•  Mg+.H20  +  C02  AH(0  K)  =  -64  kJ  mof1  (6) 

Mg02+  +  H20  -»•  Mg+.H20  +  02  AH(0  K)  =  -4 1  kJ  mof1  (7) 

Mg+.C02  +  02  -»•  Mg02+  +  C02  AH(0  K)  =  -23  kJ  mof1  (8) 

Mg+.N2  +  02  ->  Mg02+  +  N2  AH(0  K)  =  -60  kJ  mof1  (9) 

The  Mg+.X  cluster  ions  may  also  undergo  addition  reactions,  in  which  a  second  ligand  binds  to 
the  complex  in  the  presence  of  a  third  body,  M: 

Mg02+  +  02  -»•  Mg02  .02  AH(0  K)  =  -7 1  kJ  mof1  ( 1 0) 

Mg+.C02  +  C02  Mg+.(C02)2  AH(0  K)  =  -42  kJ  mol1  (11) 

Mg+.H20  +  H20  ->  Mg+.(H20)2  AH(0  K)  =  -105  kJ  mol1  (12) 

Mg02+.  +  N2  ->  Mg02+.N2  AH(0  K)  =  -96  kJ  mof1  (13) 

Mg02+  +  C02  ->  Mg02+.C02  A 7/(0  K)  =  -144  kJ  mof'  (14) 

Mg+.H20  +  02  -►  Mg02+.H20  AH(0  K)  =  -173  kJ  mol1  (15) 

It  has  been  suggested  that  formation  of  these  larger  cluster  ions  could  be  the  first  step  in  forming 
ice  nuclei  for  noctilucent  clouds  in  the  mesosphere,35  although  this  requires  the  clusters  to  grow 
very  rapidly  so  that  they  are  not  completely  destroyed  when  DR  occurs. 

Above  95  km,  where  the  chemistry  of  Es  layers  is  of  particular  interest,  MgO+  formed  in  reaction 
1  can  react  with  O3  via  two  exothermic  channels: 

MgO+  +  03  — ►  Mg+  +  202  AH(0  K)  =  -178  kJ  mol1  (16a) 

—>  Mg02+  +  02  AH(0  K)  =  -268  kJ  mof1  (16b) 

Atomic  O  is  a  major  constituent  of  the  lower  thermosphere,  and  the  sequence  of  reactions 

Mg02+  +  O  -»•  MgO+  +  02  AH(0  K)  =  -126  kJ  mol1 

MgO+  +  O  -»•  Mg+  +  02  AH(0  K)  =  -277  kJ  mof1 
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(17) 

(18) 


reduces  these  oxide  ions  back  to  Mg+.  If  reactions  17  and  18  compete  effectively  with  DR,  this 
will  limit  the  efficiency  with  which  Mg+  is  converted  to  Mg,  thereby  increasing  the  [Mg+]  /  [Mg] 
ratio  and  the  lifetime  of  Mg+  in  an  Es  layer. 

It  should  be  mentioned  here  that  the  reaction 

MgO+  +  H20  -»•  MgOH+  +  OH  (19) 

O  1 

is  fast.  However,  above  80  km  the  ratio  [O]  /  [H20]  is  so  large  (>100)  that  reaction  19  is 
uncompetitive  with  reaction  17. 

Reactions  1  and  3  have  been  studied  previously  at  300  K  using  the  flowing  afterglow  technique, 
but  the  rate  coefficients  obtained  were  suggested  to  have  an  uncertainty  exceeding  a  factor  of 
2. 30,31  Reactions  2-5  have  been  studied  at  294  K  using  a  selected-ion  flow  tube  (SIFT) 
apparatus;  however,  in  all  cases  no  reaction  was  observed.36  In  this  project  reactions  1-18  were 
studied  using  the  pulsed  laser  photolysis/laser  induced  fluorescence  (PLP-LIF)  and  fast  flow 
tube/mass  spectrometer  (FT-MS)  techniques,  where  appropriate.  We  also  studied  the  reaction 

Mg+  +  N20  -»•  MgO+  +  N2  Atf(0  K)  =  -55  kJ  moF1  (20a) 

-►  Mg+.N20  Atf(0  K)  =  -71  kJ  mol'1  (20b) 

to  investigate  the  surprising  finding36  that  Mg+  is  unreactive  towards  N20,  in  contrast  to  the  rapid 
reactions  of  other  metallic  ions  such  as  Ca+  37  and  Fe+  38  with  this  oxidant. 

In  terms  of  the  neutral  chemistry  of  magnesium  in  the  MLT,  Figure  1  shows  that  atomic  Mg  is 
oxidised  rapidly  by  O3  to  MgO,39  which  then  reacts  with  02,  O3,  H20  or  C02.40  A  key  question 
is:  how  stable  is  Mg03,  or  other  species  that  form  from  it,  and  do  they  provide  a  permanent  sink 
for  magnesium  below  90  km?  To  answer  this  question,  the  reactions  of  MgO,  Mg02,  Mg03  and 
MgC03  with  atomic  O  were  studied: 


MgO  +  O  — > 

Mg  +  O2 

AH(0  K)  = 

-235  kJ  mol 

(21) 

Mg02  +  O  — 

*  MgO  +  O2 

AH(0  K)  = 

-172  kJ  mol'1 

(22) 

Mg03  +  O  — 

»  Mg02  +  O2 

AH(0  K)  = 

-1 15  kJ  mol'1 

(23) 

MgC03  +  O 

— »  Mg02  +  CO2 

AH(0  K)  = 

-1 17  kJ  mol'1 

(24) 

The  rate  coefficients  of  these  reactions  were  measured  using  a  fast  flow  tube  with  detection  of 
Mg  and  MgO  by  laser  induced  fluorescence.  In  order  to  establish  the  technique  the  reaction 

MgO  +  CO  — »  Mg  +  C02  A/7(0  K)  =  -272  kJ  mol'1 

was  also  studied,  where  CO  can  be  regarded  as  a  surrogate  for  atomic  O. 
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Experimental  Work 


Pulsed  laser  photolysis/laser  induced  fluorescence  (PLP-LIF) 


Figure  2.  Schematic  diagram  of  the  pulsed  laser 
photolysis  /  laser  induced  fluorescence  apparatus:  fu 
flow  of  reactant  and  H e;f2,  flow  of  He;  L  =  lens  (f = 
45cm);  PMT  =  photomultiplier;  MC=  monochromator. 


The  pulsed  laser  photolysis/laser 
induced  fluorescence  (PLP-LIF) 
apparatus  in  Figure  2  has  been 
described  in  detail  by  Plane  et  a  I.41 
Mg+  was  produced  in  the  central 
chamber  of  the  stainless  steel  reactor 
by  pulsed  multi-photon  photolysis  of 
magnesium  acetyl  acetonate  (MgAcAc 
or  Mg(C5H702)2)  in  an  excess  of  the 
bath  gas  (He)  and  reactants  (CO2,  H20, 
N2,  O3  or  O2).  Powdered  MgAcAc  was 
placed  into  a  tantalum  boat,  located  in 
a  heat  pipe  connected  to  the  central 
chamber,  and  heated  to  between  433 
and  463  K,  maintained  to  within  ±  5  K 
during  a  particular  experiment.  The 
resulting  MgAcAc  vapour  was 
entrained  in  the  bath  gas  and  carried  to 
the  central  chamber  where  it  then 
mixed  with  a  larger  flow  of  bath  gas 
and  reactant  mixture. 


The  photolysis  of  the  MgAcAc  vapour 
at  193  nm  by  an  ArF  excimer  laser  (pulse  energy  40-110  mJ,  pulse  width  25  ns,  pulse  rate  5  Hz) 
produced  Mg+  ions  and  initiated  the  reaction.  The  Mg+  ions  were  probed  at  279.6  nm  (Mg+(32P- 
32S))  using  a  frequency  doubled  Nd-YAG  dye  laser  (Sirah,  Model  CBR-G-30;  laser  dye 
Rhodamine  590)  frequency  doubled  using  a  BBO  crystal  (typical  pulse  energy  =  400  pJ).  The 
excimer  and  dye  laser  beams  were  aligned  collinearly  to  pass  through  the  centre  of  the  chamber, 
with  a  dichroic  mirror  used  to  protect  the  dye  laser  from  the  excimer  beam.  The  LIF  signal  was 
measured  by  a  photomultiplier  tube  after  passing  through  an  interference  filter  centred  at  280  nm 
(fwhm  =15  nm).  The  LIF  signal  was  collected  by  a  gated  integrator,  delayed  for  30  ns  after  the 
dye  laser  pulse;  the  dye  laser  was  triggered  at  a  time  delay  after  the  excimer  pulse  that  was 
scanned  under  computer  control. 


7 


For  reaction  1,  the  O3  concentration  was  measured  using  UV  absorption  spectroscopy  between 
250  and  300  nm.  The  broad-band  radiation  from  a  deuterium  lamp  (Oriel,  Model  60000)  was 
passed  through  a  1  m  pathlength  absorption  cell,  located  downstream  from  the  central  chamber, 
and  focused  into  the  entrance  slit  of  a  0.5  m  grating  spectrometer  (Spex,  model  1870B)  equipped 
with  a  1200  groove  mm'1  grating  (resolution  0.12  nm  FWHM).  Absorption  spectra  were 
recorded  with  a  photodiode  array  (EG&G,  PARC  1412)  and  converted  to  optical  density,  before 
fitting  reference  absorption  cross  sections  for  O3  to  derive  their  concentrations.42  The  O3 
concentrations  were  corrected  for  wall  loss  and  temperature  difference  between  the  reaction 
chamber  and  absorption  cell. 


Fast  flow  tube-mass  spectrometer  (FT-MS) 


Turbopump 


Figure  3.  Schematic  diagram  of  the  fast  flow  tube  with  detection  by  mass  spectrometry  (ions)  or 
laser  induced  fluorescence  (neutrals). 

The  fast  flow  tube-mass  spectrometer  system  in  Figure  3  has  been  used  previously  to  study  both 
the  neutral  and  ion-molecule  reactions  of  Fe-  and  Ca-containing  species.43'46  The  stainless  steel 
flow  tube  (internal  diameter  =  37.5  mm)  consists  of  sections  of  tube,  cross-pieces  and  nipple 
sections  connected  by  conflat  flanges  sealed  with  copper  gaskets.  The  tube  has  a  total  length  of 
1285  mm  from  the  carrier  gas  entry  point  to  the  mass  spectrometer  skimmer  cone.  Mg+  ions 
were  produced  via  laser  ablation  of  a  piece  of  magnesite  (MgCO^)  using  a  532  nm  Nd:YAG  laser 
(repetition  rate  =  10  Hz,  pulse  energy  ~  lOmJ),  loosely  focused  onto  the  target  using  a  quartz 
lens  (focal  length  =150  mm).  The  ablation  target  was  mounted  on  a  rotary  feedthrough  powered 
by  a  DC  motor,  and  projected  into  the  centre  of  the  cylindrical  axis  of  the  tube  (not  shown  in 
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Figure  3).  The  target  was  rotated  so  that  a  fresh  surface  was  presented  to  each  laser  pulse  in  order 
to  maintain  a  uniform  Mg+  signal.  The  Mg+  ion  pulses  were  entrained  in  a  flow  of  He  which 
entered  upstream  of  the  ablation  target.  Reactants  entered  the  flow  tube  via  side  ports 

downstream  of  the  ablation  target.  An 
overall  gas  flow  rate  of  3000  seem  was 
used  at  pressures  of  1.2  -  1.7  torr, 
creating  flow  velocities  of  32  -  22  m  s'1. 
Hence  the  Reynolds  number  was 
always  less  that  80,  ensuring  laminar 
flow  within  the  tube. 

Mg+  and  product  molecular  ions  were 
detected  using  a  differentially  pumped 
2-stage  quadrupole  mass  spectrometer 
(VG  Quadrupoles,  Model  SXP  Elite). 
-200  o  200  400  600  800  1000  1200  The  entrance  aperture  of  the  skimmer 

Time  /  jus  cone  before  the  first  stage  of  the  mass 


3000 

2500 

2000 
=3 
CO 

+X  1500 

+  D) 

“  1000 
500 

0 

0.012  0.016  0.020  0.024  0.028 

Time  /  s 


spectrometer  was  0.1  mm  diameter, 
located  1160  mm  downstream  of  the 
ablation  target.  The  aperture  diameter 
of  the  skimmer  cone  before  the  second 
stage  (containing  the  quadrupole  and 
Channeltron  electron  multiplier)  was  2 
mm.  In  order  to  reduce  the  residence 
time  of  the  ions  between  the  skimmer 
zones,  the  first  skimmer  cone  was 
biased  by-llto-16V  and  the  second 
by  -94  V.  The  resulting  potential 
difference  accelerated  the  ions  into  the 


Figure  4.  Top  panel:  time-resolved  profile  of  the  LIF 
signal  obtained  by  pumping  the  (Mg+(32Pi)  -  Mg+  (3 1 S0)) 
transition  at  279.6  nm  and  monitoring  emission  at  the 
same  wavelength,  following  the  pulsed  photolysis  at 
193.3  nm  of  MgAcAc.  The  solid  line  is  a  fit  to  the  form 
A  exp1  A  n.  Temperature  =  190  K;  Pressure  =  4  Torr. 
Bottom  panel:  Diffusional  loss  of  Mg+  in  the  flow  tube  at 
295  K.  Left  axis:  Pulses  of  Mg+  for  a  series  of  flow  times 
at  a  constant  pressure  of  0.9  Torr.  Right  axis:  integrated 
signal  for  these  pulses  (filled  circles)  and  exponential  fit 
yielding  the  diffusional  loss  of  Mg+. 


quadrupole  and  prevented  the  formation 
of  clusters  ions  in  the  cold  expansion 
after  the  first  skimmer  cone.  A  multi¬ 
channel  scaler,  synchronized  with  the 
YAG  laser,  was  used  to  record  the  ion 
signal,  typically  from  the  accumulation 

'J 

of  10  laser  shots. 

For  the  neutral  experiments,  Mg  vapour 
was  produced  by  heating  the  metal  in  a 
crucible  at  about  800  K.  Mg  and  MgO 
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were  detected  by  laser  induced  fluorescence  at  285.2  nm  (Mg(3 1 P i  -  3 1  So))  and  499.4  nm 
(MgO(B1Z+  -  X'Z+),  respectively,  using  a  pair  of  Nd:YAG-pumped  dye  lasers. 


Results  and  Discussion 


PLP-LIF  measurements 

The  time -resolved  LIF  signals  in  the  pulsed  laser  experiments  were  of  a  single  exponential  form 
and  were  well  fitted  by  the  form  A.exp(-AA),  as  shown  in  Figure  4  (top  panel).  The  loss  of  Mg+  is 
described  by  the  pseudo  first-order  decay  coefficient  k',  since  the  concentrations  of  the  reactant 
and  bath  gas  were  kept  well  in  excess  of  Mg+.  For  reaction  with  O3  (reaction  1), 

k'  =  &DiffMg+  +  &MgAcAc  [  Mg  A  C  Ac]  +  /Cl  [O3]  (I) 

where  &DiffMg+  describes  the  diffusion  of  the  Mg+  ions  out  of  the  volume  defined  by  the 
intersection  of  the  laser  beams  within  the  field  of  view  of  the  photomultiplier  tube,  and  &MgAcAc  is 
the  rate  coefficient  for  the  reaction  between  Mg+  and  the  organometallic  precursor.  The  value  of 
(^DiffMg+  +  ^Mgx[MgAcAc])  was  determined  from  a  fit  of  the  LIF  decay  when  [O3]  =  0,  and 
ranged  from  3000  -  20000  s'1  depending  on  the  heat-pipe  temperature  for  a  particular 
experiment.  Plots  of  k'  versus  [O3],  which  should  be  linear,  showed  marked  curvature  at  large 
[O3].  This  behaviour  is  caused  by  reaction  16a  recycling  MgO+  back  to  Mg+.  This  problem  was 
overcome  by  adding  HC1,  which  removes  MgO+  via  reaction  2 1 : 

MgO+  +  HC1  MgCl+  +  OH  Atf(0  K)  =  -83  kJ  moF1  (26a) 

-►  MgOH+  +  Cl  AH(0  K)  =  -9 1  kJ  moF1  (26b) 

The  rate  coefficient  k\(,  for  the  reaction  between  MgO+  and  O3,  and  the  branching  ratio  /i6a  to 
form  Mg+  (i.e.,  k\eJk\k),  were  then  estimated  by  modelling  the  observed  dependence  of  k'  on  [O3] 
in  the  absence  of  HC1. 

For  the  recombination  reactions  between  Mg+  and  X  in  the  presence  of  He 

k'  =  A'D,nMg+  +  ^MgAcAc  [MgAc  Ac]  +  kx[X]  [He]  (II) 

the  rate  was  measured  both  as  a  function  of  [X]  and  [He],  Reaction  4  was  also  studied  at  high 
CO2  concentrations  to  obtain  the  relative  efficiency  of  CO2  compared  to  He  as  third  body. 
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FT-MS  measurements 


Figure  1  (bottom  panel)  illustrates  a  sequence  of  time-resolved  Mg+  pulses  measured  by  the  mass 
spectrometer  at  the  downstream  end  of  the  flow  tube.  For  kinetic  measurements,  the  relative  ion 
concentration  can  be  determined  either  from  the  pulse  height  or  the  integrated  area  of  each  pulse, 
which  gives  essentially  the  same  result.  As  an  example,  the  first-order  removal  rate  of  Mg+  in 
reaction  3  was  obtained  from  the  relation:46 

In 

k=  £3  [02]  [He]=  -— 


where  t  is  the  time  between  the  injection  point  of  O2  into  the  flow  tube  and  the  downstream 

skimmer  cone;  and  [Mg+]  c-  are  the  measured  ion  signals  in  the  presence  and  absence  of 

O2,  respectively,  t  is  calculated  from  the  measured  velocity  of  the  Mg+  pulse  travelling  down  the 
tube,  and  is  then  reduced  by  5%  to  correct  for  the  O2  mixing  by  diffusion  from  the  injection  point 
at  the  wall  across  the  centre  of  the  flow  tube.  Radial  diffusion  and  loss  on  the  flow  tube  walls  is 
significant  in  the  FT-MS  system.  This  was  therefore  measured  for  each  species  -  ions,  atomic  O 
-  and  included  in  a  model  of  the  flow  tube  kinetics. 


The  addition  of  a  second  ligand  X  to  Mg+.X  was  studied  by  monitoring  the  dependence  of 
Mg+.X  on  [X],  which  is  described  by  the  following  differential  equations: 


d[Mg+] 
d  t 


^diff,Mg+ 


[Mg+  ]  -  kx  [Mg+  ][X] 


(IV) 


d[M^X1  =  -*awx[Mg*.X]  +  ix[Mg*  ][X]  -  *x2[Mg*  ,X][X] 

where  kx  is  the  second-order  rate  coefficient  for  the  reaction  between  Mg+  and  X,  kXi  is  the 
second-order  rate  coefficient  for  the  reaction  between  Mg+.X  and  X,  and  k^ur  is  the  rate  of  wall 
loss  by  diffusion  and  uptake.  Equations  IV  and  V  were  solved  numerically  using  4th-order 
Runge-Kutta  integration  to  find  the  value  of  kx 2  which  best  fitted  the  experimental  data  points 
(using  a  y}  minimisation46).  The  statistical  uncertainty  in  kXi  was  estimated  from  the  standard 
deviation  of  the  difference  between  kxi  obtained  from  a  fit  to  each  individual  experimental  point, 
and  the  value  from  the  global  fit.  The  overall  systematic  error  was  then  estimated  by  using  a 
Monte  Carlo  procedure  which  varied  the  model  parameters  within  their  la  uncertainty  range  (104 
times)  and  repeating  a  global  fit  to  the  experimental  points.  The  ligand-exchange  reactions 
(Mg+.X  +  Y  — >  Mg+.Y  +  X)  and  the  neutral  reactions  were  studied  in  a  similar  way. 
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The  results  of  all  the  ion-molecule  work  are  summarised  in  Table  1,  and  the  neutral  reactions  are 
listed  in  Table  2. 

Table  1.  Measured  rate  coefficients  for  ion-molecule  reactions  (la  uncertainties). 


Reaction 


1  Mg+  +  03  -►  MgO+  +  02 

2  Mg+  +  N2  +  (He)  -+  Mg+.N2 

3  Mg+  +  02  (+  He)  -+  Mg02+ 

4  Mg+  +  C02  (+  He)  Mg+.C02 
4a  Mg+  +  C02  (+  C02)  -»•  Mg+.C02 

5  Mg+  +  H20  -»•  Mg+.H20  +  He 

6  Mg+.C02  +  H20  -»•  Mg+.H20  +  C02 

7  Mg02+  +  H20  -»•  Mg+.H20  +  02 

8  Mg+.C02  +  02  -»•  Mg+.02  +  C02 

9  Mg+.N2  +  02  — ■>  Mg+.02  +  N2 

10  Mg02+  +  02 — ►  Mg02.02 

11  Mg+.C02  +  C02  —*■  Mg+.(C02)2 

12  Mg+.H20  +  H20  -►  Mg+.(H20)2 

1 3  Mg02+  +  N2  -»■  Mg02  .N2 

14  Mg02+  +  C02  — >  Mg02+.C02 

1 5  Mg.H20  +  02  -+  Mg02+.H20 

16  MgO+  +  03  -»•  Mg+  +  202 

— >  Mg02+  +  02 

17  Mg02+  +  O— ►  MgO+  +  02 


k(295  K  unless  specified) 
cm3  molecule'1  s'1 
(1.2  ±0.2)  x  I0'y  (190-  340  K) 

(2.7  ±  0.3)  x  10’31  (7/300  K)'1’88  a  (190  K) 

(3.5  ±  0.3)  x  10’31  (7/300  k)('2'55±(U4)  a  (190-299  K) 

(7.1  ±  1.1)  x  10’30  (7/300  K)('1'86±0'03)a  (  190-403  K) 

(5.3  ±  0.7)  x  10'29  (7/300  k/'186*003^  (190  -  403  K) 

(2.7  ±  0.3)  x  10'29  (7/300  K)'2'58  a 

(5.1  ±0.9)  x  10’11 

(1.9  ±  0.6)  x  10'11 

(2.2  ±0.8)  x  10'11 

(3.5  ±  1.5)  x  10'12 

(3.3  ±0.3)  x  10'13b 

(8.9  ±3.8)  x  I0’13b 

(2.6  ±0.4)  x  I0’12b 

(1.9  ±0.4)  x  10'13b 

(7.4 ±  1.2)  x  10'12b 

(1.5  ±  0.5)  x  10’llb 

(3.0  ±  1.5)  x  10’10 

(5.5  ±  1.5)  x  10’10 

(6.5  ±  1.8)  x  10'10 
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18  MgO+  +  O— »  Mg+  +  02 

19  Mg+  +  N20  —*■  MgO+  +  N2 


(193- 373  K) 


(5.9  ±2.4)  x  10'10 

(5.8  ±  1.4)  x  10-30  (T/300  K)('L94±008)a 
a  Units:  cm6  molecule'2  s'1.  b  Association  reaction  at  a  pressure  of  1.2  Ton- 


Table  2.  Measured  rate  coefficients  for  reactions  of  neutral  Mg  species  (la  uncertainties). 


Reaction 

£(295  K) 

3  11 

cm  molecule'  s' 

21 

MgO  +  0  — > 

Mg  +  O2 

(6.2  ±  1.1)  x  lO'10 

22 

MgC>2  +  0  — 

>  MgO  +  O2 

(8.4  ±2.8)  x  10'11 

23 

Mg03  +  0  — 

> Mg02  +  O2 

<5  x  10'13a 

24 

MgC03  +  0 

— >  Mg02  +  CO2 

(6.7  ±  1.8)  x  10'12 

25 

MgO  +  CO  - 

->  Mg  +  CO2 

(1.1  ±0.9)  x  10'11 

a  This  measurement  is  complicated  by  the  probable  recombination  of  MgC>2  with  O2  to  form 
MgC>4,  so  the  upper  limit  describes  the  overall  reduction  of  Mg03  to  Mg. 


Theoretical  calculations 


In  collaboration  with  Prof.  Timothy  Wright  (University  of  Nottingham)  we  carried  out  a  set  of 
high-level  quantum  theory  calculations  on  the  Mg+.X  and  X.Mg+.Y  clusters.  B3LYP 
optimizations  were  performed  employing  6-3 1  l±G(2d,p)  basis  sets.  In  several  cases  a  number  of 
different  orientations  were  investigated  in  order  to  determine  the  geometries  of  lowest  energy; 
and  in  cases  involving  O  and  O2,  different  spin  states  were  also  considered.  In  order  to  establish 
accurate  dissociation  energeties,  up  to  RCCSD(T)  single-point  energy  calculations  were  also 
employed,  using  quadruple-^  basis  sets.  Figure  5  illustrates  the  optimised  geometries  of  some 
of  these  clusters.  These  calculations  provided  the  input  data  -  vibrational  frequencies,  rotational 
constants  and  energetics  -  required  for  Rice-Ramsperger-Kassell-Markus  (RRKM)  theory 
calculations  which  enabling  rate  coefficients  to  be  extrapolated  to  temperatures  and  pressures 
that  are  not  achievable  in  the  laboratory.  Figure  6  illustrates  RRKM  fits  to  the  measured  rate 
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coefficients  for  the  ion-molecule  association  reactions.  This  shows  that  some  of  these  reactions 
(e.g.  Mg+  +  H2O)  are  significantly  in  the  “fall-off’  region  between  third-  and  second-order 
kinetic  behaviour  even  at  a  pressure  of  12  Torr. 
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Figure  5.  Geometries  of  [X-Mg-Y]+  complexes,  optimized  using  the 
B3LYP/6-3 1  l+g(2d,p)  level  of  quantum  theory.  Mg  =  green;  N  =  blue; 
O  =  red;  C  =  black;  H  =  white. 
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kl  10'30  cm6  molecule2  s'1  k  / 10'  cm  molecule 


Figure  6.  Top  panel:  third-order  rate  coefficients  for  the  recombination  of  Mg+  with  H20, 
C02,  02  and  N2  as  a  function  of  temperature,  with  M  =  He.  The  black  lines  are  low  pressure 
limiting  rate  coefficients  obtained  from  fits  of  RRKM  theory  to  the  experimental  data.  The 
discrete  symbols  are  the  experimental  data  taken  at  12  Torr.  This  shows  that  the  measured 
rate  coefficients  for  C02,  02  and  N2  are  essentially  at  their  low  pressure  limits,  whereas  the 
H20  rate  coefficient  shows  significant  fall-off,  illustrated  by  the  RRKM  fit  at  12  Torr  (grey 
solid  line).  Bottom  panel:  third-order  rate  coefficients  for  the  recombination  of  Mg+.H20, 
Mg+.C02  and  Mg02+  ions  with  H20,  C02,  02  and/or  N2  as  a  function  of  temperature,  with  M 
=  He.  The  discrete  points  are  experimental  measurements  at  a  pressure  of  1.2  Torr.  The  grey 
lines  show  RRKM  fits  at  this  pressure,  and  the  black  lines  indicate  the  low  pressure  limiting 
rate  coefficients. 
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Altitude  /  km 


Satellite  retrievals 


SCIAMACHY  limb  observations  of  Mg  Charlotte  Whalley  and  the  PI  travelled 

1  -  10  July  2005  ,  TT  .  .  „  „  .  . , 

50  -  75°  N  t0  t‘ie  University  of  Bremen  in  May 

2009  to  work  with  Dr  Miriam 

Sinnhuber  on  SCIAMACHY  retrievals 

of  Mg  and  Mg+.  An  important  result  of 

this  work  is  that  atomic  Mg  does  not 

disappear  in  the  presence  of 

noctilucent  clouds.  These  ice  clouds 

form  in  the  summer  at  high  latitudes, 

at  an  altitude  of  about  83  km.  All  other 

meteoric  metals  -  Na,  Fe,  K  -  are 

rapidly  removed  in  the  presence  of  the 

clouds,  7-49  so  that  there  are 

insignificant  concentations  of  these 

metal  atoms  below  85  km.  However, 

as  shown  in  Figure  7,  there  is  no 

significant  difference  in  the  profile  of 

atomic  Mg  retrieved  from 

SCIAMACHY  when  ice  clouds  are 

present.  We  have  shown  using  quantum  calculations  that  the  reason  for  this  unexpected  result  is 

that  Mg  atoms  stick  much  less  strongly  to  an  ice  surface  -  the  binding  energy  is  only  20  kJ  mol"1, 


[Mg]  /  cm'3 

Figure  7.  Retrieved  average  vertical  profiles  of  atomic 
Mg  in  the  presence  and  absence  of  noctilucent  clouds 
(NLCs).The  broken  lines  indicate  2a  uncertainties. 


compared  to  more  than  30  kJ  mol" 
than  1  ms  even  at  130  K. 


for  these  other  metals,  so  that  a  Mg  atom  will  desorb  in  less 


Atmospheric  Modelling 


Mars 

Low-lying  plasma  layers  have  been  observed  sporadically  in  the  Martian  atmosphere  by  radio 
occultation  measurements  from  spacecraft  such  as  the  Mars  Express  Orbiter50  and  the  Mars 
Global  Surveyor. 51  These  layers  are  just  a  few  km  wide,  and  tend  to  occur  around  90  km.  It  has 
been  proposed  that  the  layers  consist  of  metallic  ions,  for  two  reasons:  they  occur  in  the 
aerobraking  region  of  the  planet  where  meteoroids  ablate;  and  they  resemble  sporadic  E  layers  in 
the  terrestrial  atmosphere  which  are  known  to  be  composed  principally  of  Fe+  and  Mg+  ions.  The 
rate  coefficients  measured  in  the  EOARD  project  were  used  to  solve  the  problem  of  how  metallic 
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ions  can  persist  in  a  CCL-rich 
atmosphere,  where  the  ions  should  be 
neutralized  rapidly  by  formation  of 
metal-C02  cluster  ions  followed  by 
dissociative  electron  recombination. 

A  model  of  magnesium  and  iron 
chemistry  in  the  Martian  atmosphere  was 
constructed,  which  includes  meteoric 
differential  ablation  rates  calculated  with 
the  Leeds  CABMOD  model,2  photo¬ 
ionization,  and  gas-phase  ion-molecule 
and  neutral  chemistry.  As  shown  in 
Figure  8,  the  model  predicts  that  nearly 
all  the  metallic  ions  between  70  and  1 10 
km  should  be  Mg+,  because  the  reactions 
of  MgC>2+  and  MgO+  with  atomic  O 
(reactions  17  and  18)  are  fast  enough  to 
prevent  these  molecular  ions  undergoing 
dissociative  electron  recombination 
(unlike  the  analogous  Fe  species).  There 
are  enough  Mg+  ions  to  form  sporadic 
layers  of  the  observed  plasma  density, 
and  the  layers  can  have  a  lifetime  against 
neutralization  in  excess  of  20  hours. 


Figure  8.  Top  panel:  modelled  height  profiles  of  Mg+, 
Mg  and  the  neutral  magnesium  reservoirs,  for  the 
globally-averaged  Martian  atmosphere  during  daytime. 
Bottom  panel:  height  profiles  of  the  corresponding  iron 
species.  Note  the  almost  complete  absence  of  Fe+, 
because  the  reactions  of  Fe02+  and  FeO+  with  atomic  O 
are  about  20  times  slower  than  the  analogous  reactions 
of  Mg02+  and  MgO+  (reactions  17  and  18). 


Earth 

Above  85  km  in  the  terrestrial  atmosphere,  reaction  1  (Mg+  +  O3)  is  the  dominant  removal 
process  of  Mg+,  although  the  resulting  MgO+  is  mostly  recycled  back  to  Mg+  by  reaction  1 8 
(MgO+  +  O).  Below  85  km,  reaction  2  (Mg+  +  N2)  becomes  more  rapid,  followed  by  reaction  3 
(Mg+  +  O2).  However,  the  resulting  Mg+.N2  will  switch  with  02  (reaction  9),  so  that  Mg02+  is 
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Figure  9.  Vertical  profiles  of  Mg+,  Mg  and  Mg(OH)2 
predicted  by  a  1-D  atmospheric  model  for  January  and  July, 
at  mid-latitudes  (40°N). 


predicted  to  be  the  major  molecular 
ion  below  85  km.  Eventually 
MgC>2+  will  switch  with  H2O 
(reaction  7),  react  with  atomic  O 
(reaction  17),  or  undergo  DR  with 
an  electron. 

This  chemistry  has  a  dramatic 
effect  on  the  lifetime  of  Mg+  in  an 
Es  layer.  The  electron  concentration 
of  an  Es  layer  is  related  to  the 
critical  frequency  (JES)  required  for 
radio  transmission  through  the 
layer.  Mg+  in  a  strong  layer  (fEs  = 

6  MHz,  [e]  =  5  x  105  cm'3)  should 
have  a  lifetime  of  ~  1  hour  at  90 
km,  but  only  6  minutes  at  85  km. 

In  contrast,  the  Mg+  in  a  weak  layer 
(fEs  <  2  MHz,  [e  ]  <  5  x  104  cm'3), 
should  have  a  lifetime  of  about  10 
hours  at  90  km.  In  fact,  comparison 
with  our  previous  work  on  the 
lifetimes  of  Fe+  44  and  Ca+  46  in 
sporadic  E  layers  shows  that  Mg+ 
has  a  much  longer  lifetime  below 
95  km.  This  arises  because 
&ig(MgO+  +  O  —*■  Mg+  +  O2)  is 
about  20  times  faster  than  the 
analogous  reactions  of  FeO+  44  and 
CaO+  ,46  The  longer  lifetime  of  Mg+ 
(by  roughly  an  order  of  magnitude) 
below  90  km  is  consistent  with 
recent  limb-scanning  satellite 


measurements  of  Mg  ,  which  show  significant  concentrations  of  this  ion  below  this  altitude. 
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The  rate  coefficients  measured  in  this  EOARD  project  (Table  1),  together  with  our  previous 
work  on  Mg  and  MgO  reactions,39’40  were  then  inserted  into  a  time-resolved  1-D  model  between 
65  and  110  km  region,  with  0.5  km  resolution.  The  new  model,  MagMOD,  is  similar  in  structure 
to  one  developed  for  sodium  chemistry  in  the  MLT,  including  the  source  s  of  the  minor 
constituents  (O3,  O,  NO+  etc.)  which  control  the  magnesium  chemistry.52  The  only  tuneable 
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parameter  in  the  model  is  the  meteoric  ablation  flux  of  Mg.  The  annual  average  flux  was  set  to 
8200  Mg  atom  cm"2  s'1,  and  this  was  then  varied  seasonally  according  to  the  radio  meteor  rate.53 
Figure  9  illustrates  the  model  predictions  of  Mg+,  Mg  and  Mg(OH)2  in  the  MLT,  for  January  and 
July  at  mid-latitudes.  This  shows  a  3-fold  increase  in  the  Mg+  column  abundance,  from  3.0  x  109 
in  winter  to  6.1  x  109  cm"2  in  summer.  This  is  in  very  good  agreement  with  the  GOME16  and 
SCIAMACHY  satellite  measurements.  The  peak  height  of  the  Mg  layer  is  between  90  and  100 
km,  with  peak  densities  between  2000  and  5000  cm"  ,  in  accord  with  rocket-borne  mass 
spectrometric  measurements.22,25 

The  neutral  Mg  layer,  in  contrast,  is  predicted  to  exhibit  little  seasonal  variation,  with  a  small 

9  2  9  2 

decrease  from  2. 1  x  10  cm"  in  January  to  1.5  x  10  cm"  in  July.  This  agrees  well  with  the 
satellite  record16,18  The  predicted  peak  of  the  Mg  layer  around  88  km  is  in  accord  with  Mg 
profiles  retrieved  from  SCIAMACHY  limb-scanning  measurements.17 

In  201 1,  MagMOD  will  be  incorporated  into  the  Whole  Atmosphere  Chemistry  Climate  Model 
(WACCM),  produced  by  the  National  Center  for  Atmospheric  Research  (Boulder).  This  general 
circulation  model,  which  extends  from  the  surface  to  140  km,  will  then  be  used  to  explore  the 
global  variation  of  Mg  and  Mg+. 
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